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Abstract 

The equations of fluid dynamics developed in paper I are applied to the study of the propagation 
of ultrasound waves. There is good agreement between the predicted propagation speed and 
experimental results for a wide range of Knudsen numbers. 
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I. STATEMENT OF THE EQUATIONS 



Modern theoretical studies of the influence of dissipation on the propagation of sound 
on the basis of the Navier-Stokes equations may be said to have begun with the work of 
Kirchhoff |T^. A principal aim of that and subsequent studies is to determine how the 



propagation speed and the rate of dissipation of the waves depend on their frequencies. For 
this problem, the predictions from the standard Navier-Stokes equations of fluid dynamics 
do not agree well with experiments when the periods of the sound waves become as short 
as the mean flight times of the particles of the gas, that is, when we enter the ultrasound 
regime. 

There are two directions from which to enter that regime. We can begin with a gas of 
freely steaming particles and introduce weak interactions among them. In that case, we 
may with Uhlenbeck pSf ask, "How is it possible to impose on the random motion of the 
molecules the ordered motion ... which a sound wave represents?" In the modern language 
of dynamical systems theory, this could be seen as a problem of synchronization in which we 
witness increasing numbers of particles going into cooperative motion until all are engulfed. 
On the other hand, we may start from the case of continuum mechanics and attempt to 
extend the validity of that description to the case of longer and longer mean free paths. It is 
unlikely that in either case we can successfully traverse the full range of possible conditions, 
but we may expect to encounter an interesting transition between the two regimes. 

In this paper, we examine how well the fluid dynamical description of paper I of this 



series [|Tn| extends into the domain where the particle mean free paths are comparable to the 
characteristic macroscopic length scale of the medium. In paper I, we derived an extension 
of the fluid dynamical equations that we hope may offer an improvement of this kind and, 
in the present paper, we study their linear form and the resultant dispersion relation for 
sound waves. In this first section, we restate the equations given in I before going on to the 
straightforward determination of the dispersion relation they imply for the linear theory of 
sound waves. 



The basic form of the macroscopic equations derived from kinetic theory, are M, & 113 



dtp + V-ipu)=0 (1.1) 
p(atu + u- Vu) + V-P = (1.2) 
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^pi?(aiT + u- VT)+P: Vu + V-Q = , (1.3) 

where p is mass density, u is the velocity field, T is the temperature, P is the pressure 
tensor, Q is the heat fiux vector and the colon stands for a double dot product. We have 
not included an external force. 

These equations express Newton's laws of motion for a continuum in phenomenological 
theories and they are a formal consequence of most kinetic theories. Where approaches to 
the derivation of these equations from kinetic theory may differ is in the expressions for the 
higher moments, P and Q. The derivations from kinetic theory are important since they 
provide formulas for the transport coefficients that appear in the specific expressions for 
the pressure tensor and the heat fiux. However, not all treatments of the kinetic theory 
give the same explicit formulas for P and Q, there being differences of degree and style of 
the approximations used. Of course, when the mean free path of the constituent particles 
is sufficiently short compared to all macroscopic lengths in the problem, there is no real 
disagreement, since the standard Navier-Stokes forms work well enough for most purposes. 
But when the macroscopic lengths become short and are comparable to the mean free 
paths of the particles, those standard results do not agree with experiment, as we shall 
see. Therefore we must ask whether there is a continuum approximation that may provide 
improved treatments of such problems. 

To test whether the expressions for P and Q derived in paper I from the relaxation model 
of kinetic theory |P, |l^ may fulfill this need, we here apply them to study of the propagation 



of ultrasound. In the relaxation model, the relaxation time, r, may be taken to be of order 
of the mean fiight time of particles, where the mean speed is of the order of the speed of 
sound. Then, we have 

^ = ^ ■ (1-4) 
where a is a constant that depends on the collision cross section and the gas constant and 
we have ignored a possible dependence of the particle cross-section time on particle speed. 

The results in I are based on an expansion in r, up to first order. Those expansions led 
to a pressure tensor, 

r /ninT 9 M 

I-/iE (1.5) 



P 

where 



p = RpT , (1.6) 



i? is the gas constant, 

E-'^^ + ^-h-uS--, (1.7) 

oxj oxi 3 

11 = TP (1.8) 

is the viscosity and D/Dt = dt + u ■ V. The resuh for /i, together with ( |1.4| ), imphes 
Maxwell's conclusion that viscosity does not depend on density for simple gases. For the 
heat flux, we obtained 

Q = -ryVT - r/TV Inp - ^/i^ (1.9) 

where r] = ^fj,R is the conductivity. Both (|1.5| ) and (|1.9| ) carry errors of order that are 
not indicated explicitly. 

These formulae for P and Q are not expressed explicitly in terms of the fluid fields. 
Rather, their expressions involve some of the same time derivatives of these fields that 
appear in the fluid equations. Here is the central difference between our results and those 
obtained in the Chapman- Enskog approach. In the latter, partial derivatives with respect to 
time are eliminated by the use of lower order results. Though we do not use that elimination 
procedure in deriving the closure relations, we can nevertheless readily recover the Navier- 
Stokes results from ours, when r ^ 0, as we described in I. However, even though both 
theories formally have first-order accuracy in r, the results from them are significantly 
different at Knudsen numbers of order unity, as we shall see in what follows. 



II. THE LINEAR THEORY 



We consider the evolution of perturbations on a uniform medium and define pertur- 
bation variables (p, 6 and w through the relations 

p = Po(l + 0), T = ro(l + ^) and p = po{l + w) , (2.1) 

where po? and po are the constant background values of the thermodynamic fields. The 
perturbation quantities all have small amplitudes with, for example, |0| << 1. From the 



linearization of ( |1.6|) we obtain 

w = (t) + e. (2.2) 

We further assume that there is no background motion so that u is small and needs no 
subscripts. 



In analogy with (|2.1|) we write for the pressure tensor, the hnearized form 



P = Pol + T 



(2.3) 



where 



T = wl - r + -V • u I - rE 



(2.4) 



where r is evaluated for the state variables of the background medium and E is small. 
Similarly, since there is no zeroth order heat flux, we get for the linearized heat flux. 



Q = -r]oToV{e + w) - -i^odtn. 



(2.5) 



The quantities /iq = rpo and rjQ are the viscosity and conductivity evaluated in terms of the 
state variables for the constant background medium. 
The linearization of ( |1 . Ij ) is 

(9i0 + V ■ u = . (2.6) 
A compact form of the linearized ( |1.2| ) is 



Po dtu + poV ■ T = 0. 



(2.7) 



When we take the divergence of (PT7|) and use (p^.6|), find that 



pod^^ - Po VV : T = 



(2.8) 



From (O) we find that 



V ■ T = - rdtVe - -rV(V ■ u) - rV • E 

o 



We see from (1/7) that 



V-E = V2u+-V(V-u). 



Then, with the help of ( p.6|) , we find 



(2.9) 



(2.10) 



V ■ T = Vtu - raV( 



rV'u. 



(2.11) 



On using ( p.6|) again, we find that 



VV : T = V^ccJ - rdtV^ie - 20) , 



(2.12) 
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which we may introduce into (pl8|). Next we define the Laplacian speed of sound, a, and the 
kinematic viscosity, z/q, as in 

a = — and z/q = — . [2.16) 
Po Po 

We then obtain the dissipative wave equation 

{d'} - a^V^ - 2uodtV^) + {uodt - a^) V^^ = . (2.14) 
To complete this discussion, it is useful to introduce the thermal diffusivity 

'^0 = ^, (2.15) 
where Cp = ^R. Thus, r/o^o = fpo^o- This is used in the linearized heat equation, 

^Podte+poV-u + V-Q = 0, (2.16) 

where we may write 

V ■ Q = -^Po«:oV'(</. + 26) + ^-po^odU ■ (2-17) 
Since z/q = a^r, we then find 

rd, - H^^v^^ 9 + ( '-r^d^ - ^-rd, - '-'^vA = (2.18) 



where 



a = ^ (2.19) 



is the Prandtl number of the undisturbed medium. 

Finally, to further simplify the appearance of these formulae, we let the unit of time be 
T and the unit of length be ar. Then our linearized equations for sound waves are 

{d^ - - 2dt\/^) (l)+idt- 1) V'e = , (2.20) 

5^2 _ _ 20^^ + (^sdt - ^V'^ 9 = . (2.21) 

For comparison we note that the analogous linear equations for the Navier-Stokes case (with 
zero bulk viscosity) are these: 

- - ^dtV^^ - V^^ = , (2.22) 
- 29(0 + (sdt - ^VM 9 = . (2.23) 
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III. THE DISPERSION RELATIONS 



We may seek solutions to the linear equations (|2.2CI|) - (|2.21|) in which and 9 vary like 



exp {ikx — st). Since the mean free path is the unit of length, the wave number k, which is 
nondimensional, is effectively the Knudsen number for this problem. The dispersion relation 
is 

(3 + 5A;2),3 - - 1^ + 5 (^h' + 1^ kh - ^fc' = . (3.1) 
For comparison, we report that the dispersion relation for the Navier-Stokes equations is 

^ ■ ' 1 ' ,2 „2 , r- / 4 2 , 1 \ , 2 „ 5 4 



3s' - (^^ + 4J k's' + 5 \^-^k' + Ij k's -^k^ = 0. (3.2) 
A. Frequencies 

To get a feeling for what these results mean, we look at free modes for which k is real. 
Then we set s = iu + a where a and u are also real. When we introduce this into ( p.l|) 
we find that there is a (thermal) mode with u = and a pair of (acoustic) modes whose 
frequencies satisfy 

2 2 2(10-0) 2 , 5((T+5P) ,^2 /oox 

which gives the frequencies of sound waves. As we may confirm, a is of order unity for large 
k and it grows in proportion to k"^ for small k. Hence, for both very large and very small 
k, the last term on the right of (|3.3| ) is the largest one on that side. So we may write the 
uniform approximation 

For small k, this gives the phase speed uo/k = ± a/5/3, which is the usual speed of sound 
for an adiabatic sound wave, as is to be expected for very long wave lengths. For large fc, 
we obtain the phase speed uj/k = -^^hjo. 

For the N-S equations with zero bulk viscosity there is the same number of modes: a 
thermal mode with zero frequency and sound waves with 

= 3c,2 _ ^ ^^^2^ + ^(1^ + 1)A;2 . (3.5) 
3 cr 3 3(7 

As expected, the two sets of equations agree in the limit of very small fc, where the N-S 



equations return the phase speed ±a/5/3. But for large /c, the differences between the two 



theories become quahtative. With the Navier-Stokes equations, we find that at large k, 
instead of reaching a finite limit, the phase speed is proportional to k for large k. As we 
shall see when we look at the experimental results, the N-S prediction is qualitatively wrong; 
the phase speed goes to a finite value at large k. 



B. Damping Rates 



The equation for the damping rate is 
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a 



3(3 + 5P)cu2 - 5(1 + -A;2)A;2 



«+(H_i)A;V--A;^ = 0. (3.6) 
a o 



For the thermal mode, for which = 0, we find the damping rates a = k'^/a for small k 
and a = 1/5 for large k. Thus, there is very little damping for long waves while short waves 
are damped on the collisional time scale. Moreover, on examination of these two limits, we 
see that they each emerge from the balance of the same two terms in ( |3.6| ). Hence we may 
write the approximate formula 

a = — (3.7) 



a + 5P 

as a reasonable approximation to the damping rate for all k, in the thermal mode. 

Similarly, in the case of sound waves, we see that a is also the result of the balance 
between the same two terms in (|3.6| ) in the limits of for large and small k. Hence, we find 
that for sound waves, to good approximation, the damping rate is given by 

(10 - a)k^uj^ - 5k^ 



a 



(3.8) 



3a(3 + 5P)u;2 _5p(^ + 5^2) ' 

where uj is given in ( |3.4| ). For long wave lengths, the damping is again slight since it goes 
to zero like [(7 — cx) / {6a)]k'^ . For large k, we obtain the finite limit a = (5 — cr)/(5cr). 
For the Navier-Stokes equations, the damping rate is given by 



3a^ - {- + A)k^a' + 



a 



3(7 



a + (- + 4)A;V--fc^ = . (3.9) 
a a 



The damping rate for the acoustic modes also goes like k"^ for small k for both the thermal 
mode and the acoustic modes. For the acoustic modes, the damping rate is a = [(2cr + 
1) / {3a)]k'^, so that we get the same wave number dependence, but with a different coefficient 
than is obtained from our equations in the small k limit. However, for increasing k, the N-S 
damping rates grow like k"^ for sound waves, which is in disagreement with experiment Pj. 



8 



IV. COMPARISON WITH EXPERIMENT 



Though the study of free modes in the previous section is intuitively clear, it does not 
directly represent the way experiments on sound propagation are usually carried out. In the 
experiments, it is more typical that one drives the fluid at a real, fixed frequency and then 
studies the propagation of waves in space. The forcing may be accomplished by vibrating the 
end wall of a tube containing gas at a fixed (real) frequency uj and observing the propagation 
down the tube. To model this procedure in full detail would involve a careful treatment of 
the forcing procedure, which usually requires attention to boundary conditions. However, 
in this first reconnaissance of the way our equations describe sound waves, we shall adopt 



a standard theoretical practice [0] and simply fix the wave frequency, a;, in the dispersion 
relation and compute the resulting fc, which will typically be complex. Thus, in (|3.1|) we let 
s = iuo and we find that the equation for k becomes 



- (1 - biu) + 



a 



hiu}'^ — hiuj — (— — 
a 



k^ + 3iuj-^ 



0. 



(4.1) 



We may similarly obtain such an equation for the N-S case, ( |3.2| ). In order to emphasize 
the results for large Knudsen number we plot the results in the manner used, for example, 
by Cercignani That is, we introduce the quantity 

5P 



(4.2) 



where i^' is a normalized inverse propagation speed. The factor 5/3 is included so that the 
phase speed is nondimensionalized on the Laplacian (or adiabatic) speed of sound, rather 
than the Newtonian (or isothermal) speed of sound as above. Then we find that 



;i - 5zcu) + 



a 



hi \ uj — 



5crz 



0. 



(4.3) 



To see how this representation contains the results for free modes, we note that, in the 
limit — > 0, (^]3D reduces to i^^ = 1. That is, for low frequencies, the usual adiabatic sound 
speed is recovered. In the (more interesting) opposite limit, lo oo, we obtain = cr/3 
for the propagative modes. Thus we see that, in the limit of forcing at high frequency, sound 



waves propagate with phase speeds ±a/o"/3, independently of frequency. The data shown in 
the accompanying figure (Fig. |l|) confirm this independence of frequency (or wave number) 
of the speed of propagation of ultrasound. The observed nondimensional phase speed is 
0.47. 
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The Prandtl number found from the relaxation model of kinetic theory, either by the 
methods of Chapman and Enskog or those described in paper I, is unity. With this value, 
we obtain 0.51 for the limiting phase speed, so this represents a small quantitative error. 
However, the value of a found in kinetic theory depends on the atomic model used, that 
is, on the nature of the collision term. Though the relaxation model gives the explicit 
value unity for a, the value found with the traditional Boltzmann equation for hard spheres 
is 2/3. This difference has nothing to do with the approximation method (our procedure 
gives (T = 2/3 when applied to the Boltzmann equation) but is a consequence of the nature 
of the form of the atomic interactions that is adopted. We therefore follow a common 
practice put the empirical Prandtl number into the theoretical results when comparing with 
experiments. Since the experimental data we shall refer to are for noble gases whose values 
of a are 0.6 or 0.7 we shall here adopt the value a = 2/3 suggested by the Boltzmann 
equation. When we use that value of the Prandtl number in evaluating the phase speed, we 
obtain K = v^/3 = 0.471. Even without this adjustment, the results for the propagation of 
ultrasound are good, but we would propose to anyone thinking of using our equations from 
this first-order development from the relaxation equation to introduce this phenomenological 
improvement of the theory. 

In the accompanying figure (Fig. |I|), we show the variation of 3?e as a function of 1/u; 
from a number of sources. The experimental values ([|I^, 0) are indicated as individual 
points (the diamonds) and they appear to be tending toward a nonzero constant value at 
high frequency. This is qualitatively in accord with our results, here shown as a solid line 
for the case of cr = 2/3, and it is in stark disagreement with the prediction from the Navier- 
Stokes equations (long dashes with double dots), which predict that 3?e K goes to zero like 
l/uj. Since our limiting value for 9fte K was found to be a/o"/3, the remarkable agreement 
of our results with experiment owes something to our using the experimental value of 2/3 
for a for limiting value K ^ 0.47. Nevertheless, even without this choice, the results 
would be adequate and comparable to those shown for the moment method (ll^) with 
16, 215 moments (short dashes). Other theoretical studies of ultrasound are based on direct 
solution of the Boltzmann equation |^, |23| and we show the results of Sirovich and Thurber 
20| obtained in this way (medium dashes), for which the Prandtl number automatically has 



the value 2/3. 



10 




10-2 10-1 10° 10' 

Inverse of Knudsen Number (cox)-^ 

FIG. 1: The inversed phase velocity as a function of inversed frequency 
V. CONCLUSION 

In the study of the thermal damping of sound waves by electromagnetic radiation , 
one finds that, for thermal times much less than the acoustic period, sound propagates at 
the isothermal speed of sound with negligible dissipation. In the opposite limit of long 
thermal times, there is also little dissipation, but propagation is at the adiabatic speed of 
sound. The experiments, and the solution of the Boltzmann equation show similar behavior 
when the relevant parameter is the ratio of the collisional relaxation time to the acoustic 
period. Our equations, as well as those of the moment method (with tens of thousands of 
moments), reproduce this behavior but the Navier-Stokes equations do not. Moreover, when 
the Prandtl number is chosen to be that of the experimental gas, the quantitative agreement 
becomes very good. In the next installation of this series, we shall compute the profile of a 
stationary shock wave. As we shall see, the agreement with the experiments is good in that 
case too. 
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